arXiv:1507.05090v2 [astro-ph.CO] 8 Dec 2015 


Mon. Not. R. Astron. Soc. 000 , (2015) Printed 10 December 2015 (MN style file v2.2) 


Weak lensing by galaxy troughs 
in DES Science Verification data 

D. Gmen^’^*, O. Friedrich^’^, A. Amara^, D. Bacon'^, C. Bonnett^, W. Hartley^, 

B. Jain®, M. Jarvis®, T. Kacprzak^, E. Krause^, A. Mana^’^, E. Rozo^, 

E. S. Ryrf’^ S. Seitz^’^ E. Sheldon^®, M. A. Troxe0\ V. Vikrami^, T. M. 

C. Abbott^^, F. B. Abdalla^'^dS^ g_ Allam^®, R. Armstrong^^, M. Banerji^^ds^ 

A. H. Bauer^®, M. R. Becker^^i^ A. Benoit-Levy^'^, G. M. Bernstein®, 

R. A. Bernstein^^, E. Bertin^^d4^ g L Bridle^^, D. Brooks^'^, E. Buckley-Geer^®, 

D. L. Burke^’®, D. GapozzE, A. Garnero Rosell^®d6^ ]yj Garrasco Kind^^^s^ 

J. Garretero®do^ M. Grocce^®, G. E. Gunha^, G. B. D’Andrea'^^s^ L. N. da 

Gosta2®d6^ D. L. DePoy^®, S. Desadb32^ H. T. DiehE®, J. P. Dietrich3b32^ 

P. DoeO^, T. F. Eifler®’^^, A. Fausti Neto^®, E. Fernandez®, B. Flaugher^®, 

P. Fosalba^®, J. Frieman^®’^"^, D. W. Gerdes^®, R. A. Gruendl^^^s^ q Gutierrez^®, 

K. Honscheid^®’^^, D. J. James^^, K. Kuehn^^, N. Kuropatkin^®, O. Lahav^'^, 

T. S. LP, M. Lima2®d9, M. A. G. Maia2®’26, M. March®, P. Martini^®^®, 

P. Melchior3®’37, G. J. Miller^®^!, R. Miquel®d2, J. J. Mohr2’3h32^ ^ Nord^®, 

R. Ogando^®d®, A. A. Plazas^^, K. Reil®, A. K. Romer^^, A. Roodman^’®, M. Sako®, 

E. Sanchez^'^, V. Scarpine^®, M. Schubnell^®, I. Sevilla-Noarbe^^d4^ (g_ gmith^^, 
M. Soares-Santos^®, F. Sobreira^®^®, E. Suchyta^®’^^, M. E. G. Swanson^®, G. Tarle^®, 
J. Thaler^®, D. Thomas^d®, A. R. Walker^^^ R. H. Wechsler^sdi^ j. Wellerh2.32^ 
Y. Zhang^®, and J. Zuntz^^ 

Author affiliations are listed at the end of this paper 
★ E-mail: dgruen@usm.uni-muenchen.de 


ABSTRACT 

We measure the weak lensing shear around galaxy troughs, i.e. the radial alignment of 
background galaxies relative to underdensities in projections of the foreground galaxy 
held over a wide range of redshift in Science Verihcation data from the Dark Energy 
Survey. Our detection of the shear signal is highly signihcant (10-15cr for the smallest 
angular scales) for troughs with the redshift range 2 ; S [0.2, 0.5] of the projected galaxy 
held and angular diameters of 10 arcmin ... 1°. These measurements probe the connec¬ 
tion between the galaxy, matter density, and convergence helds. By assuming galaxies 
are biased tracers of the matter density with Poissonian noise, we hnd agreement of 
our measurements with predictions in a hducial A cold dark matter model. The predic¬ 
tion for the lensing signal on large trough scales is virtually independent of the details 
of the underlying model for the connection of galaxies and matter. Our comparison 
of the shear around troughs with that around cylinders with large galaxy counts is 
consistent with a symmetry between galaxy and matter over- and underdensities. In 
addition, we measure the two-point angular correlation of troughs with galaxies which, 
in contrast to the lensing signal, is sensitive to galaxy bias on all scales. The lensing 
signal of troughs and their clustering with galaxies is therefore a promising probe of 
the statistical properties of matter underdensities and their connection to the galaxy 
held. 
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1 INTRODUCTION 


The measurement of weak gravitational leasing probes 
matter inhomogeneities in the Universe by means 
of the differential deflection they induce on the 
light of background sources. Most leasing analyses 
are driven by the signatures of matter overdensi- 


ties, e.g. the gra v itational shear of galaxies (e.g. 


Brainerd et al.1 19961: Hoekstra et al.1 l2004l: ISheldon et al. 

2OO4I: iMandelbaum et al.l 20061: Ivan Uitert et al. 

2OIII: Brimioulle et al.1 l2013l: 

Velander et al.1 20l4 

Glamoitt et al.l in DreoaratiorJl and clusters of Ecalax- 

ies (e.e. Tvson et al.1 1990l: 

iMarrone et al.l 20091: 

Sheldon et al. 20091: iMarrone et al 

2012: Hoekstra et al. 

2OI2I: von der Linden et al.l 12014 

: iGruen et al.1 I20l4 

Umetsu et al. 20141') or the spatial correlation of shear 


due to intervening large-scale structure, cosmic shear (e.g . 


Wittman et al. 

2000l:lFu et al.ll2008 

:ISchrabback et al.ll201(]i 

Kilbineer et al 

2OI3I: Becker et al.l 

201,5l'l. 


Probes of underdense structures are complementary to 
this. Differences between dark energy and modified grav¬ 
ity (MG) models for cosmic acceler ation might be more 
easily differentiable in cosmic vo ids l|ClamDitt et al.ll2013l : 
iLam et ^ l2014l : ICai et ahl l2014h . The reason for this is 
that the screening of the hypothetica l fifth force of MG 
(|Vainshteinl 1 19721 : iKhourv fc WeltmanI 1200^ 1. required to 
meet Solar system constraints, is absent in these low-density 
environments. MG therefore entails that negative density 
perturbations should grow more rapidly than predicted by 
General Relativity, with eff ects on the dens ity profile in and 
around such structures (cf. ICai et al.llioi^ b 


Using voids detected in Sloan D igital Sky Survey 
(SPSS) spectroscopic galax y catalogues JSutter et al.ll20ia . 
I2OI4I : iLeclerca et al.112OI5I I , first measurements of the ra- 
dial alignment of background galaxies hav e been made 
dMelchior et al.l[2014l : Oampitt fc Jain|[2014l b Future spec¬ 
troscopic surveys will yield lensing measurements of void 
matter profiles wi th moderate signal-to-noise ratio (SNR; 

I Krause et al ]l2013lb Comb ined with predictions for void pro¬ 
files ( Hamaus et al ]|2014al ia. these will provide unique tests 
of gravity. 

On large enough scales, we expect a symmetry between 
the excess and deficit of matter relative to the mean. De¬ 
spite this fact, highly significant lensing measurements have 
thus far only been performed on matter overdensities. In this 
work, we follow the new approach of measuring the proper¬ 
ties of underdense regions in projections of the galaxy den¬ 
sity field over wide ranges in the radial coordinate, i.e. in 
redshift. Due to the wide redshift range (e.g. z € [0.2, 0.5]) 
used for the projection, these can be straightforwardly iden¬ 
tified in galaxy catalogues with photometric redshift (photo- 
z) estimates. Because the selection from the projected field 
avoids underdensities which are randomly aligned with mas¬ 
sive structures in front or behind them along their lines of 
sight, the SNR of the lensing signal is comparatively high. 
The approach is related to and motivated by the measure¬ 
ment of shear aroun d galaxies in und e rdens e projected envi¬ 
ronments, for which iBrimioulle et al.l (|201.ll 'l have previously 
detected significa nt radial aligment (cf. their fig. 25; see also 
iGillis et al.l201 j for a detection of radial shear around galax¬ 
ies in underdense 3D environments). 

We make these measurements using Science Verifi¬ 


cation (SV) data from the Dark Ene r CT Su rv ey (DES, 
The Dark Energy Survey CollaborationI I2OO5I : iFlaughen 
2OO5II . The data was taken after th e commission i ng of 
the Dark Ene rgy C amera (DECam; IPePov et al.1 120081 : 
iFlaugher et al.1 l2015h on the 4m Blanco telescope at the 
Cerro Tololo Inter-American Observatory (CTIO) in Chile 
to ensure the data quality necessary for DES. We make use 
of a contiguous area of 139 deg^ for wh ich SV imaging data 
at lensing quality is available (cf. also IVikram et al. 20151 : 
I Jarvis et al.ll20f^ . 

The structure of this paper is as follows. In Section 2, 
we describe the data used. Our theoretical modelling of the 
trough signal is introduced in Section 3. Section 4 presents 
our measurements of the shear signal around troughs and 
their angular two-point correlation with galaxies. We sum¬ 
marize and give an outlook to future work in Section 5. For 
all theory calculations, we use a fiducial fiat A cold dark mat¬ 
ter (ACDM) cosmology (flm = 0.3, fib = 0.044, us = 0.79, 
Us = 0.96). 


2 DATA 

We select trough positions and measure their lensing signal 
with galaxy catalogues from the SV phase of DES. In this 
section, we briefly describe the catalogues used. 


2.1 Galaxy catalogue 

The DES SVAl r ed-sequence Matched Filter Galaxy Gat- 
alog (redMaGiC; IRozo et al.1 l2015lj is a photometrically 
selected luminous red galaxy (LRG) sample chosen to 
have precise and accurate photometric redshifts. red- 
MaGiG makes use of the red sequenc e model computed 
from the redMaPPer cluste r catalogue iiRvkoff et al.ll2014l : 
iRvkoff et al.l l in preparatiorJl . This model of the red sequence 
as a function of magnitude and redshift is used to compute 
the best-fitting photo-2 for all galaxies regardless of SED 
type, as well as the goodness-of-fit. At any given red¬ 
shift, all galaxies fainter than a minimum luminosity thresh¬ 
old Lmin are rejected. In addition, redMaGiG applies a x^ 
cut such that x^ < Xmax, where the maximum Xmax('Z) is 
chosen to ensure that the resulting galaxy sample has a 
nearly constant comoving space density n. In this work, we 
use the high density sample, such that h = 10“®/i®Mpc“® 
and I/min(2:) = 0.5I/*(2) in the assumed cosmology. This 
space density is roughly four ti mes that of t he SP SS BOSS 
CMASS sample. As detailed in iRozo et al.l ll2015ll . the red¬ 
MaGiG photo-2s are nearly unbiased, with a scatter of 
(JzIG + z) ~ 1.7% and a 4a outlier rate of about 1.7%. 
By virtue of this, redshift errors are negligible for the se¬ 
lection of underdensities in the galaxy field when the latter 
is projected over the wide redshift ranges we use (see Sec¬ 
tion [2I2Il- 


2.2 Trough selection 

The troughs are selected as centres of cylindrical regions (or, 
more accurately, of conical frustums) of low galaxy density 
as follows. Let the galaxy catalogue (cf. Section r2.1ll be given 
with entries Xi,Zi,Li for the angular position, redshift and 
luminosity of galaxy i, 1 ^ i ^ n. 
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Define a function M4(2, L) that assigns a weight to each 
galaxy based on its redshift and luminosity. Furthermore, 
define Wt{0) to weight points by their projected angular 
separation 9 from the centre. From these we define a pro¬ 
jected, weighted and smoothed version of the galaxy density 
field 


i=l 


In this study, we use a simple weighting corresponding 
to a hard cut in luminosity, redshift and radius, i.e. 


WW,L) = 
Wt{9) = 


1, 

L ^ 0.5L* 

0, 

otherwise 

1, 

9 ^ Ot 

0, 

otherwise 


A 0.2 s: 2 s: 0.5 


( 2 ) 


with some trough radius Ot, which we vary between 5 and 
30 arcmin. The redshift range of 2 € [0.2, 0.5] is motivated 
by a trade-off between having a large trough volume (and 
hence a large signal) and having sufficient background galax¬ 
ies with useful shear estimates to measure the effect, al¬ 
though we test other settings in Section [4.1.11 Note that a 
luminosity-independent weighting scheme of LRGs, as we 
apply it here, is no t far from optim al for the reconstruction 
of the matter field llCai et al.ll201ll i. 

Probed on a finely spaced grid of sky positions, for 
which we use an Nside = 4096 healpix dOorski et al.llioosl i 
map with 0.86 arcmin pixel spacing, there is a distribution 
of weighted galaxy counts as measured by G. We select the 
set of trough positions T as the points below the 20th per¬ 
centile G 20 of that distribution, i.e. T = {x : G{x) ^ G2o}EI 
It is important to note that we are not selecting individual 
minima of the galaxy density field. Instead, troughs are over¬ 
lapping within and between samples of fixed trough radius. 
A map of trough positions for = 10 arcmin is shown in 
Fig. [T] Note how the troughs visibly trace underdensities in 
the convergence field. 

In practice, G needs to be corrected for the effect of 
masking due to survey boundaries or bright stars, for ex¬ 
ample. A homogeneous masking fraction over all cylinders 
simply decreases the tracer density, which is accounted for 
by our model automatically. When masking fractions are not 
homogeneous but vary among the cylinders, the situation is 
more complicated. For the present analysis, we make the 
approximation of excluding all cylinders where more than 
/ = 20% of the area inside the trough radius is masked in 
the galaxy catalogue and assume that all remaining cylin¬ 
ders have equal masking fractions when selecting troughs 
and modelling the signal. At the level of statistical precision 
achieved here, this simplification is not expected to cause a 
significant difference. 

The mean surface density of redMaGiC galaxies in the 
useable area is approximately fi — 0.055 arcmin“^, corre¬ 
sponding to a mean count of approximately 4, 17, 69 and 
155 galaxies in cylinders of radius 5, 10, 20 and 30 arcmin. 


^ We have also tried stricter (i.e. lower percentile) and looser 
(i.e. higher percentile) thresholds, which yield a higher (lower) 
amplitude of the signal with larger (smaller) uncertainties. Also 
see Section 14.1.21 for the selection of high-dens\ty cylinders. 



Figure 1. Positions of 10 arcmin troughs (black circles, to scale, 
randomly selected sample of 1500 out of the fiill0,000 trough 
positions probed) overlaid on to lensing convergence map (red: 
p ositive, blue: negative ). Convergence was estimated as described 
in IVikram et ahl |201^ with a weighted lensing source catalogue 
fSection l2.3ll and a 7 arcmin Gaussian smoothing. 


At the lower 20th percentile, selected troughs have mean 
counts of 1, 9, 44 and 108 galaxies, respectively. 


2.3 Lensing source catalogue 

For the background sources, we use a shape catalogue 
measured with NGMIX0 We apply the cuts, w eighting and 
respon sivity correction as recommended in I Jarvis et al.1 
(l2015l) . where also the shape measurement and testing of 
catalogues is described in detail. In order to prevent confir¬ 
mation bias, shear estimates in the catalogue were blinded 
wit h an unknown fact or until the analysis had been finalized 
(cf. Liar vis et al ] l2015l. their section 7.5). 

We use the two highest redshift bins defined in 
iBecker et al.l ll2015l . cf. their fig. 3) by means o f photo- 2 : prob¬ 
ability density estim ates made with SkyNet llGraff fc Ferod 
l2013l: Ib onnettl[2015l ). a method t hat performed well in an 
extensive set of tests on SV data llSanchez et al.llioidl ). The 
mean redshift of the lower (higher) redshift bin is 2 : « 0.6 
{z « 0.9). We use the appropriately weighted SkyNet 
stacked p[z) estimate for predicting the lensing signal (cf. 
Section [ 3 ]). To maximize the SNR for our non-tomographic 
measurements, we weight the signal measured in both bins 
as 1:2 to approximately accommodate the ratio of effective 
inverse critical surface mass density. The resulting p{z) of 
the samples used are shown in Fig.[2l The tests performed in 
iBonnett et al.1 (l2015l ) indicate that errors in photo-z estima- 


^ cf. https;//github. com/esheldon/ngmix 
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Figure 2. Redshift distributions of the low-redshift leasing source 
sample (blue, short dashed line), high-redshift sample (red, long 
dashed line) and the combined fiducial source sample (black, solid 
line). All distributions are weighted by inverse shape noise of 
the sources in analogy to the shear signal and normalized for 
J p(z) dz = 1. Redshift distribution of all redMaGiC galaxies is 
shown as the magenta dotted line (solid for the fiducial trough 
redshift range z G 0.2,..., 0.5 indicated by the dotted vertical 
lines). 


tion are not a dominant systematic error for the prediction 
of the shear power spectrum that we use them for. 

We note that the consistency of different shape 
measurement and photo-z methods with the catalogues 
used in this analysis has been investigated in detail 
and confirmed within the systema tic requirements on 


the present data in several works | 

Sanchez et al.l 

2014 

Jarvis et al.l 20151: Bonnett et al. 2015 

: Becker et al.l 

2015 

The Dark Enerev Survev Gollaboration et al.ll2015l). 


3 THEORY 


Structures in the Universe can be described by an underly¬ 
ing field, the matter density as a function of position and 
time. Matter density itself is not an observable. Its prop¬ 
erties can, however, be recovered by a number of observ¬ 
able fields, such as the three-dimensional or projected galaxy 
density (a sparse, biased tracer) or the convergence field (a 
weighted, projected version of it). In this section we describe 
our modelling of the shear and galaxy correlation signal of 
troughs by the interrelation of these observables and the 
underlying matter density. 

We assume that the three-dimensional galaxy field can 
be described as a deterministic, biased tracer of the matter. 
This means that the 3D contrast <5 of matter density p, 


<5 = 


P-{P) 

(P) 


(3) 


and the equivalent quantity defined for the galaxy field are 


proportional at any position. Their ratio defines the bias b, 
which depends on the galaxy population. 

Lensing convergence k is related to <5 via the 
projection integral over com oving distance x (cf-, e-g- 
iBartelmann fc Schneideill200ll ). 

poo 

= / dX <5(X^',X) , (4) 

Jo 

where 

o 

(5) 

with 

poo / 

G(x) = / dx' nBource(x') ^ ■ (6) 

Here, nsource(x) is the distribution of source galaxies of the 
lensing measurement. 

An overdensity of convergence inside a circular aper¬ 
ture relative to its edge results in a tangential alignment of 
background galaxy shapes. Correspondingly, an underden- 
sity causes radial alignm ent. Both cases are described by 
(cf. ISchneider et al.l[2006l . p. 279f) 

-ytie) = (k)(< 6) - k,{6) . (7) 

Here, 7t(S) is the tangential component of gravitational 
shear averaged over the edge of a circle of radius 9, k{6) 
is the equivalent average of the convergence and {«)(< 9) 
is the mean convergence inside the circle. For the case of 
|7t| "C 1 , |k| -C 1 , tangential components of gravitational 
shear and reduced shear gt are approximately equal and ob¬ 
servable as the mean tangential alignment of background 
galaxy ellipticity. The cross component of shear, 7x, rotated 
by 45° relative to the tangential direction, is expected to be 
zero when taking the average over the full circle for a single 
thin lens or over an ensemble of thick lenses. 

In order to connect these fields and model the trough 
signal, we make these three assumptions: 

• We apply the lLimberj (Il954l) approximation to compute 
the angular power spectrum of the projected matter density 
contrast 6^(0) within the redshift range of the redMaGiC 
galaxies used for the trough selection (0.2 ^ x ^ 0.5). The 
same approximation is also used to compute the cross power 
spectrum between 5-e and the convergence field k{0) relative 
to the background galaxy redshift distribution. 

• We assume that &(0) and k.{0) follow a Gaussian dis- 
tribntion - at least when they are averaged over the trough 
radius or over the annuli in which we measure the shear. 

• We assume that the redMaGiG galaxies are placed on 
to (5 e via a biased Poisson process. 

In this section, we describe how these assumptions trans¬ 
late to a prediction for the expected shear signal and galaxy 
density around troughs. 


3.1 Projected matter density and galaxy counts 

Let the volume density of redMaGiC galaxies as a function 
of comoving distance x be given by niens(x)- The projected 
galaxy contrast is proportional to a weighted projection of 
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matter contrast (5 e. In a flat universe, the latter is calculated 
as 

^ ^ ^ = [ dx niens(x) , (8) 

^ •ixo 

where S{x0,x) is the 3D matter contrast at the point 
{X0, x) 0^1 til® backward light cone. For a galaxy sample 
with constant comoving density, such as the redMaGiC cat¬ 
alogue, this is a simple volume weighting of matter density, 
Wiens (x) dy/[dfldx], and fe is the common projected 
matter density contrast. 

If we average (5 e over circles of angular radius 0 t, we 
arrive at the new random field dx, 

i / d'e' ■ (9) 

™T J|6i-0'|<eT 

In this we have made the approximation of a flat sky, valid 
for 6t <C 1. 

If the galaxies are placed on to fe {&) via a biased Pois¬ 
son process, then the discrete probability P of finding N 
galaxies within Ot given the value of dx is 

P(A^|5x) = ^ (Njl + Mx])"^ exp(-iV[l-bMx]) . (10) 

For ^x < we assume P{N > 0|(lx) = 0 and P{N = 
0|5x) = 1 (see also Appendix A). We have used the bias b 
and mean galaxy count N within 6t- For our model predic¬ 
tions shown later, we fix the bias at a fiducial value of fe = 1.6 
or vary it between 1.4,..., 1.8 to show the dependence on 
bias in the relevant range. Note that we neglect the mod- 
era te redshift depen dence of the bias of redMaGiC galaxies 
(cf. IRozo et al.lboisl f. which is a good approximation for the 
limited redshift range used here. 

We identify troughs as circles in the sky with low galaxy 
count N. Given any N, the expected value of (5x is 


((5x|Ai) = / d^x (Ix p((Ix|Al) . 


( 11 ) 


Bayes’ theorem tells us that 


p{St\N) = 


P(iV|<5x)p(5x) 

P(fV) 


( 12 ) 


When the variance Var((Ix) := <C 1, we can approximate 
p((5x) by a Gaussian distribution, i.e. 


p{St) 



(13) 


In Appendix B, we derive how can be calculated 
from the power spectrum of fe, which in the lLimberj (Il954| j 
approximation is related to the 3D matter po wer spectrum 
Ps by (cf., e.g. iBartelmann fc Schneideill20^ . Eqn. 2.84) 


Ge(£) = ^ dx Ps(^^,x^ ■ (14) 

For our calculatio ns, we use the no n-linear matter power 
spectrum Ps from ISmith et all ll2003l ). 

The normalization factor in equation 1121 P{N), gives 
the probability of finding N galaxies inside a cylinder of ra¬ 
dius 0x. Fig. [3] shows the observed distribution of galaxy 
counts in cylinders. There is reasonable agreement to our 
prediction from a Gaussian random field 5x on to which 



1 10 100 


N 


Figure 3. Distribution of redMaGiC galaxy counts in zx S 
[0.2,0.5] for various trough radii Sx- We show measurements 
(black, solid line) and model predictions for Gaussian matter 
fields on to which galaxies are placed as biased tracers with 
b = 1.6 and Poissonian noise (red, long-dashed lines). Model 
predictions for a lognormal matter density are also shown (blue, 
short-dashed lines). Normalization of P{N) is matched to make 
max(P®‘^“®®'‘^“) = 1 for each trough radius. 


galaxies are placed by means of equation [TO] (dashed, red 
lines). The distribution of counts inside 30 arcmin radii ap¬ 
pears more peaked than the model, potentially explained by 
the fact that there is too little area to have enough uncor¬ 
related troughs at this scale. 

As an alternative prescription, we also calculate the ex¬ 
pected P{N) for the case of a lognormally distributed mat¬ 
ter density. This is done by replacing the Gaussian p(5x) of 
equation [T3] by a shifted lognormal distribution of the same 
variance with a minimum value of 5x = —1. The result (dot¬ 
ted, blue lines in Fig. 0 resembles the low-density tail of the 
galaxy count distribution more closely. 


3.2 Convergence and shear profile 


Our goal is to compute the mean value of k at a distance 6 
from the trough centre. For this, we define a set of annuli 
i = 1,..., n around the trough for which d £ [6i, 6i+i). Let 
Ki be the average of k in annulus Ai, 


If both 5x and Ki have Gaussian distributions with zero 
mean, then the expectation value of Ki for a fixed value of 
5x is given by 


= s) 


Gov(<5x,R0 


a 


2 

X 


S , 


(16) 


where the covariance Gov((5x, Ki) can be computed in terms 
of the cross power spectrum Ck,s{£) of E and k (see 
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Appendix B). Note that the Gaussian approximation for 
the matter contrast and convergence becomes accurate re¬ 
gardless of the pointwise p{S) since all random fields are 
smoothed over annuli or circles and a large redshift range. 

The expectation value of Ki when the trough contains 
N galaxies is finally given by 


{K,\N) = / ds = s> p(5t = s\N) 


Cov{ST,Ki) 


ds s p(St = sjN) 


Cov(St,K,) 


(<5T|iV> . 


(17) 


If we select as troughs all cylinders with N ^ N^r 
galaxies, the mean Ki around them will be 


{Ki\ ^ Amax) = 


Cov(^t,ao 

4 nN) 


(18) 


The tangential shear signal around troughs is given by equa¬ 
tion [7] using equation [T5] to calculate the mean convergence 
in each annulus. 


3.3 Trough-galaxy angular correlation 


The trough-galaxy angular correlation function can be mod¬ 
elled in a very similar way. First define annuli At around 
the trough that correspond to the bins in which w{6) is 
measured. The mean density contrast Wi in each annulus is 
given by (cf. equation [TS]) 


Under the assumptions of Gaussianity, the expectation value 
of 6i for a fixed value of St is given by (cf. equation 1161) 


= s) 


C0Y{ST,Wi) 


(7 


2 

T 


S . 


( 20 ) 


In analogy to equation [TS] the mean density contrast in an¬ 
nulus Ai around the trough is given by 


(Wil < Amax) = 


Cov{ST,Wi) EnToPW {St\N) 

'^T EnTo PiN) 


■ ( 21 ) 


The average number of galaxies in an annulus i outside 
the trough radius is given by 


A- 

(A, I ^ An,ax)out = + H■W^\ S: Aniax)] 

= Ai[l + b{Wi\ < Amax)] , 


( 22 ) 


where the mean galaxy count Ai in annulus i is obtained by 
rescaling the average galaxy count inside one trough radius 
to the area Ai of the annulus, 

Ai = A^ . (23) 

The profile of galaxy counts around a trough is then given 
by 

{wnM ^ Amax)out = ^ _ 1 

= b{Wi\ ^ Amax) . (24) 


The situation is more complicated for annuli inside the 
trough radius Ot- Here, the Poisson noise of the different 
bins is correlated. This is because the sum of the galaxy 
counts in the different bins has to meet the requirement by 
which we selected the troughs. 

Without full treatment of the covariances, we make a 
prediction for the galaxy number counts inside the trough 
that (i) matches the mean galaxy counts predicted by the 
full model and (ii) matches the projected matter contrast 
profile with a given bias. To this end, we simply replace A in 
equation[22]by the predicted mean number of galaxies inside 
selected troughs, e.g. when demanding that At ^ Amax. 
This number is given by 


At 


T.nZ 7 PjN) A 
EvT PiN) 


(25) 


The mean number of galaxies found in an annulus inside the 
trough is then given by 

(Ail =5 Amax)in = At [1 + 6(Wi| Sj Amax)] (26) 
J\.T 


and the profile of galaxy counts inside the trough radius is 
given by 


(tUv.i] ^ Amax)in 


(Ai] ^ Amax)in _ ^ 

Ai 

:^(l+&(M;i] S; Amax)]-1 . 


(27) 


4 MEASUREMENT 

In the following section, we correlate trough positions with 
the shear signal of background galaxies fSection 14.11) . Ad¬ 
ditionally, we measure the projected number density profile 
of redMaGiC galaxies in the same redshift and luminosity 
range used for the trough selection, i.e. the angular two-point 
cross-correlation of troughs and galaxies f Section 14.21) . 


4.1 Shear signal 

We measure the mean shear of background galaxies around 
troughs, selected as described in Section [2.21 To correct for 
potential additive shear systematic errors, we subtract the 
tangential shear measured around random points. Since the 
masked region depends on the respective trough radius, the 
random shears for each 0 t differ slightly. Fig. U shows mea¬ 
sured tangential and cross shears. 

Per-mille radial alignment of background galaxies at 
and beyond the trough radius is detected with high signif¬ 
icance in all bins (cf. Section [4.1.31 Table 1). Gross shears 
are consistent to the expected null signal within the uncer¬ 
tainties (cf. also the reduced Xx in Table 1). The model 
proposed in Section[2]is a good fit to the data in all bins (cf. 
Section and reduced xEa in Table 1). 

4-1.1 Tomography 

By splitting either the source sample or using smaller red¬ 
shift ranges for selecting the troughs, it is possible to probe 
the redshift evolution of the trough leasing signal. We per¬ 
form both measurements in the following. 
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Figure 4. Weak lensing signal of galaxy troughs of =5, 10, 20 and 30 arcmin radius (top left to bottom right). Shown is the 
tangential shear signal (blue) around points from the lower 20th percentile in galaxy counts in cylinders of z = 0.2,..., 0.5. Lines show 
model predictions (cf. Section 0 for our fiducial cosmology and, for illustration of the bias dependence, a bias of 6 = 1.4,1.6,1.8 (light 
to dark blue, dotted, dot-dashed and dashed lines). Cross-shear is shown with grey cross symbols, to be interpreted with error bars of 
similar size. Tangential shear around random points, subtracted from the trough measurement, is shown with black open symbols. 


• For source tomography, we divide the source galaxy 
sample into two redshift bins (cf. Section 12.311 . Note that 
since troughs are thick lenses, the change in source redshift 
causes more than a simple change in amplitude. The dif¬ 
ferential weighting as a function of lens redshift inside the 
2 = 0.2 ... 0.5 cylinder also influences the shape of the shear 
profile. Due to the nearly power-law matter two-point cor¬ 
relation at all redshifts, however, the latter effect is small. 
The left-hand panel of Fig. [5] shows the source-tomographic 
signal. We note that the agreement of the measurement with 


the model in both bins is additional evidence for the appro- 
pri ateness of the p(z)s as estimated for our source samples 
(cf. iBonnett et al.l[20^ ). 

• For trough redshift tomography, we split the trough 
redshift range into two approximately equal-volume slices 
z = 0.2 ... 0.4 and z = 0.4 ... 0.5. When using these smaller 
redshift ranges for the trough selection, two effects reduce 
the SNR: (1) due to the lower galaxy count, Poissonian 
noise weakens the correlation of trough positions with mat¬ 
ter underdensity; and (2) uncorrelated, overdense large-scale 
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Figure 5. Source tomography (left) and trough redshift range tomography (right) of trough lensing signal. Left-hand panel shows 
tangential shear signal for sources in the lower (blue, short dashed) and higher (red, long dashed) of our redshift bins (cf. Section 12.311 . 
Right-hand panel shows signal of troughs selected by the galaxy count in zjp = 0.2,..., 0.4 (blue, short dashed) and z = 0.4,..., 0.5 (red, 
long dashed). Model predictions (cf. Section[3ll are shown for a bias of 6 = 1.6 and grey/black points indicate gx for both measurements. 



6 [arcmin] 



6 [arcmin] 


Figure 6. Tangential shear signal for troughs, i.e. centres of cylinders below the 20th (measurement and model in blue), and overdense 
cylinders above the 80th percentile (red) in galaxy count. We plot model predictions for a bias of b = 1.6, with solid (dashed) lines 
assuming a Gaussian and dotted (dashed-dotted) lines a lognormal distribution of the matter contrast around troughs (overdense 
cylinders). Grey/black points indicate gx for both measurements 


structure along the line-of-siglit outside the trough redshift 
range causes additional variance in the lensing signal. Shear 
measurements are shown in the right-hand panel of Fig. [S] 
The signal is reduced as expected, but the measurement is 
still highly significant and consistent with the model in both 


cases (see Table 1 for details on significance and goodness 
of fit). 
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4-1.2 Galaxy density percentiles 

All measurements presented above use troughs selected to 
be below the lower 20th percentile of galaxy couirts. Mea¬ 
surements with larger limiting percentiles (e.g. the 30th per¬ 
centile) give results of similar significance but smaller am¬ 
plitude. 

It is particularly interesting, however, to study the sym¬ 
metry of matter in the overdense and underdense tails of 
the galaxy field. For dense enough tracers and large enough 
scales, the expectation is that all involved fields are approx¬ 
imated well by a Gaussian distribution. This should lead to 
symmetric shear signals at the same upper and lower per¬ 
centiles. On smaller scales, the galaxy couirts (if only due to 
Poisson noise) and the matter density and convergence field 
(since |(5| -C 1 is no longer true and non-linear evolution 
boosts high-density fluctuations) deviate from a Gaussian 
distribution and we expect some degree of asymmetry be¬ 
tween the low- and high-density signal. 

The measurement for both the lower and upper 20th 
percentile is shown in Fig. and is in agreement with these 
expectations. At small trough radii, there appears to be a 
significant asymmetry, with the overdense regions showing 
a larger shear signal than anticipated from our model or the 
measurement of underdensities. For larger cylinders, such an 
effect is not detected. A lognormal model of the matter con¬ 
trast (dashed lines) makes virtually no difference for larger 
trough radii. For smaller trough radii, the shears around 
high-density cylinders are predicted to be somewhat larger, 
yet not sufficiently so to fit the data well. We hypothesize 
that the discrepancy between high- and low-density cylin¬ 
ders can rather be explained by an environment dependence 
of the bias of the redMaGiG tracer galaxies: because the 
mean bias of galaxies in overdense regions is larger than 
in underdense regions, the shear around small, high-density 
cylinders gets boosted relative to the signal around the low- 
density troughs (cf. the bias dependence of the model pre¬ 
diction in Fig. n. 


- 0.4 

- 0.2 

- 0.0 

r 

■ - 0.4 

Figure 7. Correlation matrix Rij of shear around 10 arcmin 
troughs measured in the logarithmic angular bins of Fig. as 
estimated from 100 jackknife regions. 
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4-1-3 Significance 

For estimating uncertainties, we use a set of Nj = 100 jack¬ 
knife resamplings. In order to ensure that these are approxi¬ 
mately equally populated with troughs, we choose them with 
a K-means algorithm^ on the catalogue of 5 arcmin trough 
positions. The delete-one jackknife yields a covariance 

Cov(/l,/ 2 ) = - {/l»(/2.- - (/2)) (28) 

■’ i=l 

for two quantities fi , /2 estimated from the data excluding 
region i (f^i) or averaging over all ((/) = N~^ /-'*)• 

our case, we estimate the covariance matrix C of tangential 
shear measurements (or, in Section 14.2[ angular two-point 
correlation measurements) in our set of angular bins. 

Fig. El shows the correlation coefficients Rij = 
Cov{gl, gi)/ \/Var(gj)Var((;t) estimated for our fiducial 
10 arcmin trough measurement. At intermediate and large 
radii, neighbouring bins are highly positively correlated, 

® https; //github. coni/esheldon/kmeaiis_radec/ 


Table 1. Metrics of significance of detection of shear around 
troughs of radius selected from the galaxy field in the given 
redshift range tx at the percentile threshold "P for sources in the 
indicated Zs bins. We list the SNR of shear at the trough radius, 
5t((?T)/o'gj(6i.p), and of the optimally weighted linear combination 
of shears, F/ot. See description in Section 14.1.31 for details. The 
remaining columns show the reduced of the residuals of model 
and measurement (cf. Section 14.311 and of cross-shears. 

which is even more the case for the larger troughs. The neg¬ 
ative correlation of the innermost bins is a generic feature 
that appears in all trough sizes probed and is connected to 
the opposite sign of the first two data points of the lower 
panels of Fig. 21 Both this and the off-diagonal negative cor¬ 
relations at large radii are also seen in less noisy versions of 
the covariance determined from simulations (cf. Friedrich et 
al., in preparation). 

We ensure the significances defined below are stable un¬ 
der a change of binning scheme and jackknife regions by cal¬ 
culating them with 15 instead of 25 radial bins for which we 
estimate the covariance using 50 rather than 100 jackknife 
patches, which yields consistent results. 
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Different measures of detection significance can be de¬ 
fined as follows: 

(i) SNR of shear. A simple measure is the tangential shear 
at the first angular bin outside the trough radius 9t in units 
of its standard deviation according to the jackknife estimate, 

9t/crg--\gt{GT)\/(^gt(err)- _ 

(iij For optimal signal-to-noise (cf., e.g. lOruen et al.l 
I 2 OIII. their Eqn. 11), we define a linear combination of tan¬ 
gential shear measurements. The weights of the linear com¬ 
bination are chosen as W oc C where we use the 

model prediction for our fiducial bias of fe = 1.6. The 

SNR of r = W ■ gt is given as T/ur = V/\/w'^ C W. 

(iii) We do not list a significance based on the of the 
null hypothesis here for two reasons: (I) since the signal is 
consistent with zero on a range of small and large scales, as is 
also expected from the model, the p-values of the measured 
strongly depend on which bins are used and (II) yields 
an uncertain estimate of SNR due to the variance Var(x^) = 
2 nbins for ribins bins. 

We list these metrics for various trough selections in Table[T] 
For the most conservative metric, gt/ug, we find a detection 
significance of lOcr for the smallest troughs. The optimal 
linear combination of observables yields even higher signifi¬ 
cances. Our detection of radial shear around underdensities 
on small scales of 6t = 5,10 arcmin is of considerably higher 
sig nificance than that of the most recent void le nsing stud¬ 
ies llMelchior et al.l[2014l : iCTampitt fc Jain|[2014ll . On larger 
scales, significance decreases but is still comparable. 

4.2 Trough-galaxy angular correlation 

The lensing signal around troughs, studied in the previous 
sections, measures a weighted, projected version of the mat¬ 
ter density field (cf. equation |T|. Galaxies themselves also 
trace the matter field, yet with a different weighting. We 
have approximated the connection of galaxies to the mat¬ 
ter density, so far, as being constant comoving density, de¬ 
terministic, biased tracers. Measurements of the two-point 
correlation of trough positions with galaxies are complemen¬ 
tary to trough lensing, sensitive to both the properties of the 
matter field and the details of the connection of galaxies and 
matter. 

We measure the angular two-point correlation between 
trough positions and redMaGiC galaxies in the same red- 
shift range of z = 0.2 ... 0.5 and limited to the same survey 
subarea of 139 deg^ also used for the lensing analysis. Un¬ 
certainties are again estimated by estimating the two-point 
correlation in 100 jackknife resamplings and are highly cor¬ 
related between bins, as is common in clustering analyses. 

Fig. [H] shows results for the fiducial trough parameters, 
i.e. the trough catalogues also used in Fig. 2] The low galaxy 
count level inside the trough, due in part to the selection of 
regions of low matter density and to Poisson noise, steeply 
rises at the trough radius outside of which there is no Poisson 
contribution. Physical, smaller underdensities in the galaxy 
field are observed out to large radii. Section 13.31 discusses 
our modelling of the signal. Although only at moderate sig¬ 
nificance, there are indications of an increase of bias with 
trough radius, related to either a gene ral scale or density de - 
pendence of bias or assembly bias (e.g. lWechsler et al.l2006li . 


4.3 Comparison to theory 

We briefly compare our measurements to the model put for¬ 
ward in Section [31 

Our measurements of tangential shear around under- 
dense troughs are consistent with the predictions at all scales 
and source and trough redshift configurations tested here. 
The reduced Xmod of fh® residual of the data with respect 
to the b = 1.6 model are listed in Tableland consistent 
with noise. It is worth noting that the model is a good fit 
essentially without any free parameter. The only exception 
to this is the mild dependence on the assumed galaxy bias 
for the smallest scale = 5 arcmin considered here. This 
can be understood as an effect of the importance of Poisson 
noise relative to true variations in the matter density field as 
traced by the galaxies. On large scales where Poisson noise 
is subdominant, galaxies are dense tracers of the smoothed 
matter field. Independent of the details of the galaxy place¬ 
ment model, i.e. as long as galaxy and matter density are 
somewhat positively correlated, the selection of some per¬ 
centile in galaxy count then yields an essentially equivalent 
selection in matter density. 

The model also consistently predicts our measurements 
for the two-point correlation of troughs and galaxies. The 
estimation of goodness of fit is strongly affected by the cor¬ 
relation of errors over wide ranges of scales. On the larger 
scales of 6t = 20 and 30 arcmin and when taking into ac¬ 
count the full covariance, the model is a good fit to the data 
for the full range of bias b = 1.6,..., 1.8 probed. For = 5 
and 10 arcmin, the large linear bias model with b = 1.8 is 
excluded at 4a significance (reduced x^ = 4.3 and 2.1, re¬ 
spectively, for the 25 data points), while the other models 
are good fits. Similarly, for the shear around overdense cylin¬ 
ders, the only deviation from the prediction is found for the 
smallest scale cylinders, where the measured shear around 
overdense cylinders is somewhat larger than predicted. Both 
observations indicate that at these scales, the linear bias of 
our galaxy placement model and/or the assumed Gaussian- 
ity of the matter field may be not completely valid. 

The combination of these measurements therefore is 
sensitive to both the details of how galaxies follow matter in 
low- and high-density environments (e.g. a scale or density 
dependence of galaxy bias) and on cosmological parameters. 
These aspects will be studied in more detail in Friedrich et 
al. (in preparation). 


5 CONCLUSIONS 

We have presented the measurement of per-mille level radial 
gravitational shear of background galaxies and negative two- 
point correlation of foreground galaxies around underdense 
cylinders (troughs) in the foreground galaxy field in DES 
data from the SV period. 

• Our detection of radial shear around these projected un¬ 
derdense regions (cf. Section r4.il) is highly significant (above 
lOcr; cf. Section [4.1.3l) . on the smallest projected scales and 
widest projection redshift range considered. This is a much 
higher significance than has been achieved with present data 
for the shear signal around three-dimensional voids. 

• We develop a model for the shear profile (cf. Section (S)), 
based on the assumption that galaxies are biased, Poisso- 
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Figure 8. Angular two-point correlation of trough positions and redMaGiC galaxies in z = 0.2... 0.5 for the same configurations as in 
Fig. a Shown are signal (black) and model predictions, for illustration for different values of the bias (b = 1.4,1.6,1.8 from light to dark 
blue, dotted, dot-dashed and dashed lines, cf. Section 0. 


nian tracers of the Gaussian matter density field. The model 
predicts the lensing measurements consistently within the 
present level of uncertainty. It is interesting to note that 
on sufficiently large scales, the prediction is virtually inde¬ 
pendent of the details of the galaxy placement model, yet 
sensitive to cosmological parameters (cf. Friedrich et ah, in 
preparation). 

• Tomographic measurements that split the source sam¬ 
ple or the redshift range used for the selection of troughs 
show consistent results. We note that the significance of 
radial shears strongly decreases for smaller trough redshift 
ranges, due to both the increased noise in galaxy counts and 


the variation of uncorrelated (overdense) structures along 
the line of sight in front or behind the trough cylinders. 

• We measure the shear signals around underdense and 
overdense cylinders in galaxy count at the same percentile 
thresholds (cf. Section r4.1.2l) . On small scales we find indica¬ 
tions for some deviation from our simple model predictions 
for the high-density regions. On large scales, however, we re¬ 
cover the expected symmetry between radial and tangential 
shear for both cases. 

• In addition to the shear signal, we measure and model 
the two-point correlation of galaxies from our tracer popula¬ 
tion around trough positions icf. 14.21) . While consistent with 
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our prediction on sufficiently large scales, this probe is more 
sensitive to the details of how galaxies trace the matter and 
therefore complementary to the shear signal. 

The statistical power of these measurements will 
strongly increase as larger data sets become available. We 
note, in particular, that the final survey area of DES will be 
« 30 times larger at comparable or even better data quality, 
allowing very precise measurements of the trough lensing 
signal. With these better statistics, trough lensing will be a 
relevant probe of cosmology, not only in the sense of con¬ 
straining parameters of a ACDM model. Also, the potential 
lack of screening mechanisms in underdense environments 
would influence the growth of negative density perturba¬ 
tions, with implications for constraining MG models with 
these measurements. 

On small scales, the details of how galaxies trace mat¬ 
ter and the intrinsic distribution of the fields involved are 
likely to play a significant role for model predictions, and 
simulations in combination with progress on modelling will 
be required. Under these prerequisites, trough lensing mea¬ 
surements are a promising tool for probing the connection of 
galaxies and matter and gravity in the underdense Universe. 
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APPENDIX A: CONDITIONAL 
PROBABILITIES 

The conditional probability density of Jt, given that N 
galaxies were found inside the radius 9t, is given by 

^ l {N[l + b5T]f 1 -4 

w m 

(29) 

where we have made the simple assumptions that galaxies 
trace matter with a constant bias b and that the variation 
of galaxy counts around the expectation value is given by 
the Poisson distribution. The normalisation constant gives 
the overall probability of finding N galaxies inside Ot, 

N = P{N). (30) 

In Appendix B, the trough variance cr^ is derived from the 
2D power spectrum of the projected matter contrast. 

Note that in order to self-consistently define the biased 
Poisson model as explained above, equation [^Hl can only be 
valid for Jt > — 1/b. Furthermore, one has to assume that 

P (N|Jt ^ -1/5) = 0 for A > 0 . (31) 

As a consequence one also has 

p(Jt|A) = 0 for A > 0 , Jt ^ —1/5 . (32) 

The case A = 0, however, is more subtle. Here one has 
P(A = 0|5 t —1/5) = 1, and hence 

p(5T|A = 0) = -p^^^for Jt^-1/5. (33) 

This also has to be considered when the probability P{N = 
0 ) is computed. 


APPENDIX B: VARIANCE AND COVARIANCE 
OF CONVERGENCE AND Jt 

Let Si, i = 1,2, be two line of sight projections of the matter 
density contrast J, i.e. 

noo 

5i[e) = j dx q^(x) b{xO,x) , (34) 

Jo 

Qi being the weights of the projectio ns (cf. 
iBartelmann & Schneiden l200ll L According to the iLimbeJ 
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lll954h approximation, the 2D cross power spectrum of 
and 82 is given by 


Cl,2(^) = 



gi (x)g2(x) 


Ps 



(35) 


Here, x is the comoving distance and a flat universe was 
assumed. Let Ai be annuli with minimal radius 0i,min and 
maximal radius Si.max- The annulus-averaged versions of 5i 
are given by 


A(6») 


/ d^e' G,[e-e') 5,{e') 

7T(6'?+l,max - ^?,max) 


(36) 


where Gi{6) is the top-hat filter corresponding to annulus 
Ai. 

5i{0) can be expanded into spherical harmonics as fol¬ 
lows: 


5i(e,^) = Y,^tiYr{e,4>) ■ (37) 

i,m 


If Si is a homogeneous and isotropic ra ndom field th en the 
coefficients satisfy the equation (cf. IPeeblesll 19931') 


) = 5u'5mm'Ci,e , (38) 

where Ci,^ is the 2D power spectrum of Si. 

Since the expectation values (Si) vanish, the covariance 
{D 1 D 2 ) can be computed as 


{D 1 D 2 ) = / dDidD 2 Gi(6li)G2(6»2)(<5i(6li)&(6l2)) 


= EE/ dDidD2 Gi(6»i)G2(6»2)(a7ar5" > = 

Yr{ei)Y,7^'{62) 

= I Gi(0i)yr(0i)x 

/? _ 


i,m 


J dD2 G2(6l2)Tr"(02) 

i Li 2 ,i , 


(39) 


wher e in the last step we used the relation (see e.g. I Peebles! 

I 1993 II 


/r = / dD f{0)Yf"'ie) . 


(40) 


The annuli and circles we will use as filters are isotropic, i.e. 
Gi{d,(f>) = Gi{6). Hence all coefficients G7r vanish except 
for G^i =: Gi,r . These are given by 

Gi,i = j dn G,(6I)F/(6I) 

In 

dcf) s\n{6)Gi{6)Pi(cos{9)) 
d6» sin(6i)Pf(cos(6i)) 

min 

dx Pe{x) . (41) 


= Nt i d0 

/•Si 


/o 


2txNi 

2 Td\fi 


i,min 
cos Bi 


Ai 


Here, Pe are the Legendre polynomials, Ai is the area of the 
annulufl and A/^ is a normalization factor given by 

The covariance then reads 

(D 1 D 2 ) = '^ChGi,,G 2 ,^ . (43) 

£ 

Note that, using the equation 

one can simplify the Gi/ to 

= (2l^ [^^+^(*) - W • (45) 

The covariance of Ki and can then be computed by 
setting 5i = Ki and S 2 = St. The variance is found by 
setting both (5i and S 2 to (5 t. 
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